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ABSTRACT 

We measure the contribution of galaxy disks to the overall gravitational potential of 30 nearly face- 
on intermediate-to-late-type spirals from the DiskMass Survey. The central vertical velocity dispersion 
of the disk stars (c^'fjLo) is related to the maximum rotation speed (Knax) as cr^'ljLo ~ 0.26V^nax, 
consistent with previous measurements for edge-on disk galaxies and a mean stellar velocity ellipsoid 
axial ratio a = a^l ^^r = 0.6. For reasonable values of disk oblateness, this relation implies these galaxy 
disks are submaximal. We find disks in our sample contribute only 15% to 30% of the dynamical mass 
within 2.2 disk scale-lengths (Jib)-, with percentages increasing systematically with luminosity, rotation 
speed and redder color. These trends indicate the mass ratio of disk-to-total matter remains at or 
below 50% at 2.2 h^. even for the most extreme, fast-rotating disks (V^ax > 300 km s^^), of the reddest 
rest-frame, face-on color {B — K ^ A mag), and highest luminosity {Mk < —26.5 mag). Therefore, 
spiral disks in general should be submaximal. Our results imply that the stellar mass-to-light ratio 
and hence the accounting of baryons in stars should be lowered by at least a factor of 3. 
Subject headings: galaxies: kinematics and dynamics — galaxies: stellar content — galaxies: halos — 
galaxies: spiral — galaxies: formation — galaxies: fundamental parameters 



1. INTRODUCTION 

Rotation-curve decomposition is the primary tool for 
measuring the distribution of dark matter in spiral 
galaxy halos, but the tool is blunted by uncertainties 
in the mass-to-light ratio (T) of the luminous disk and 
bulge. The 'maximum-disk hypothesis' (van Albada & 
Sancisi 1986; see also Binney & Tremaine 2008) by- 
passes this mass-decomposition ('disk-halo') degeneracy. 
A disk contributing maximally to the gravitational po- 
tential sets a lower limit on the amount of halo dark mat- 
ter in the inner regions of disk galaxies. Maximum-disk 
decompositions find the disk mass produces 85 ± 10% of 
the observed rotation velocity at 2.2 disk scale- lengths 
{hji; Sackett 1997). Unfortunately, this hypothesis re- 
mains unproven, and there have been suggestions to the 
contrary based on the lack of surface-brightness depen- 
dence in the TuUy-Fisher relation (TF; TuUy & Fisher 
1977) for a wide range of spirals (Courteau & Rix 1999; 
Courteau et al. 2003; Zwaan et al. 1995). Constraints 
on T from stellar population synthesis (SPS) models are 
poor due to long-standing uncertainties in the low-mass 
end of the stellar IMF and late phases of stellar evolu- 
tion (e.g., TP-AGB stars; Maraston 2005, Conroy et al. 
2009). 

Direct kinematic evidence for the mass contribution 
from spiral disks stems from a series of longslit spec- 
troscopic studies (van der Kruit & Freeman 1984; Bot- 
tema 1993; Kregel et al. 2005) showing that the ra- 
tio of the disk stellar velocity dispersion and the max- 
imum observed rotation speed was much lower than ex- 

^ University of Wisconsin, Department of Astronomy, 475 N. 
Charter St., Madison, WI 53706; mab@astro.wisc.edu 

^ University of Groningen, Kapteyn Astronomical Institute, 
Landleven 12, 9747 AD Groningen, Netherlands 

^ International Research Fellow, National Science Foundation 

^ NRG Herzberg Institute of Astrophysics, 5071 W Saanich 
Road, Victoria, BG V9E 2E7, Ganada 

^ National Optical Astronomical Observatories, Tucson, AZ 



pected for a maximum disk. Bottema (1993) found 
disks contribute only 63 ± 10% of the observed rotation 
speed. Most of these galaxies are viewed edge-on, how- 
ever, so that the inference of disk mass requires some 
uncertain assumptions about the stellar velocity ellip- 
soid (SVE) and line-of-sight deprojection. Recent results 
from PNe kinematics of several nearby galaxies suggest 
disk-maximality may depend on Hubble type (Herrmann 
& CiarduUo 2009). Constraints on disk maximality from 
hydrodynamical modeling of observed non-axisymmetric 
gas flows has conflicted (cf. Weiner et al. 2001 and Kranz 
et al. 2003). Analyses of one gravitional-lens system 
(Dutton et al. 2011) and one barred, resonance-ring sys- 
tem (Byrd et al. 2006) indicate the disks in these two 
galaxies are submaximal. Observational evidence sug- 
gests disks are submaximal, but the question is unset- 
tled. 

The DiskMass Survey (DMS; Bershady et al. 2010a) 
breaks the disk-halo degeneracy by obtaining indepen- 
dent measures of the total dynamical mass and dynam- 
ical disk-mass surface density from integral-field stellar 
and gas kinematics of face-on, disk-dominated galaxies. 
Following the approach of Bottema (1993) and Kregel et 
al. (2005), in this letter we derive the disk maximality 
based on the relation between the disk vertical stellar 
velocity dispersion (dz) and the disk-equatorial circular 
speed of the potential, making basic assumptions of disk 
equilibrium and a testable hypothesis about the SVE. We 
adopt ffo = 73 km s"^ Mpc"^ and Ve ga magnitudes. 

2. DATA AND MEASUREMENTS 

For 30 nearly face-on spiral galaxies covering a range 
in morphology, mass, color, surface brightness and scale 
length from our Phase-B sample (Bershady et al. 2010a), 
the stellar line-of-sight velocity dispersions ctlos and gas 
rotation curves were measured out to ~ 3/i/j. Spectro- 
scopic data in the Mgl region (498-538 nm) at resolutions 
of A/ A A ~ 7700 were collected with the PPak integral 
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field unit (Verheijen et al. 2004; Kelz et al. 2006) of 
tfie PMAS spectrograpli (Roth et al. 2005) on tlie 3.5m 
Calar Alto telescope^ Stellar velocity dispersions were 
measured down to ctlos ~ 18 km for individual fiber 
spectra via our data-censored cross-correlation software 
(DCS; Westfall et al. 2011a) using a Kl III stellar tem- 
plate. 

The maximum rotation speed of the potential (Knax) 
at 2.2/ijj is derived by fitting a simple hyperbolic-tangent 
(tanh) model (e.g., Andersen et al. 2001) to the two- 
dimensional [OIIIJA5007 velocity field, deprojected us- 
ing inclinations (itf) inferred from inverting the TuUy- 
Fisher relations from Verheijen (2001) for the if-band. 
To minimize systematic error, total magnitudes estima- 
tion matched the protocol used by Verheijen (2001), us- 
ing elliptical-aperture surface-photometry of reprocessed 
2MASS images (Skrutskie et al. 2006), as described in 
WestfaU et al. (2011b; DMS-IV). Since Knax is derived 
from the asymptotic velocity of the tanh model, system- 
atic errors arise if the rotation curve is not truly flat. In- 
spection of each galaxy's position- velocity diagram shows 
these errors to be of the same order as the formal errors 
from the fitting, in most cases. In Figure 1, below, we 
flag two exceptions: UGC 1862, which still has a rising 
rotation curve at the outermost point; and UGC 4458, 
which has a declining rotation curve after reaching Vmax 
due to a prominent bulge. Random errors on Vmax (7% 
on average) are dominated by inclination uncertainties. 

The central value of the stellar velocity dispersion of 
the disk (cloI) is found by fitting a radial exponential 
function to the ctlos measurements in a radial range that 
shows an exponential decline (typically 0.5 < R/hj^ < 4), 
e.g., excluding the bulge region, and then extrapolating 
to i? = 0. The exponential function is found to fit well in 
this radial range. Because CTz is expected to scale with the 
square root of the surface density, the kinematic scale- 
length (her) is expected to be 2hji for an exponential 
disk of constant thickness and T viewed perfectly face- 
on. In our fitting, we allow h^, to be a free parameter; 
the implications for T disk gradients are discussed in 
Martinsson (2011). For present purposes it suffices to 
state T gradients are generally small. 

The vertical component tT^"*^ of the stellar velocity el- 
lipsoid is extracted from af^^^ by assuming a constant 
shape of the SVE and irp- Following Bershady et al. 
(2010b; DMS-II) we adopt a = 0.6 ± 0.15 (25%) and 
13 = ag/aji = 0.7 ± 0.04 (5%), where a_R and ae are 
the radial and tangential components of the SVE. These 
are reasonable values given extant results for external 
galaxies summarized by van der Kruit & de Grijs (1999), 
Shapiro et al. (2003), and our own work (Westfall 2009 
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and DMS-IV). We may then write = ^^^i ^ ; 7 is a 
projection factor defined in DMS-II as a function of the 
SVE ratios and inclination. A typical value of y/j cos i is 
1.1, with a 10% uncertainty. This factor dominates the 
error on at small radii and larger inclinations, with 
the largest contribution from the uncertainty in a. Er- 
rors contributed by inclination uncertainties are lower 
but non-negligible; errors contributed by /3 computed 
from the epicycle approximation are negligible. Mea- 

^ Centro Astronomico Hispano Aleman (CAHA) at Calar Alto, 
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surement error in ulos generally dominate at larger radii 
and at lower inclinations. Details of the sample, instru- 
mental configurations, data acquisition and reduction, 
basic data products, and full mass decompositions are 
presented in Martinsson (2011). 

3. DISK MAXIMALITY 

For an oblate, self-gravitating, exponential disk of con- 
stant T it can be shown (Freeman 1970, Casertano 1983, 
Kuijken & Gilmore 1989, DMS-II) that the maximum 
rotation speed of the disk is related to the disk central 
mass surface-density, Eqi and /i^,: 



disk 
max 



x/tt G Y.0 hn 
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where Cmax ^ 0.88(1 — 0.28/(7) to excellent approximation 
for g < 4, and q = hji/hz is the oblateness parameter. 
A comparison of 14iax to V^jf^ is an accurate measure 
of disk maximality regardless of a bulge component. The 
comparison also provides a good estimate of baryon max- 
imality because V^'^lfx occurs at R/hn ~ 2.2 where most 
bulges do not significantly contribute to Knax- 

For a disk in equilibrium Edyn = cr'^./irkGhz, where k 
parameterizes the vertical density distribution with likely 
values between 1.5 to 2 (exponential to isothermal dis- 
tributions, respectively; van der Kruit 1988). Because 
we measure ctlos outside of the galaxy center, where the 
bulge contribution is minimal and it is plausible that 
the SVE is relatively constant with radius, our estimate 
of ct^'IILo is well-posed. In contrast, central measure- 
ments of (Jlos are contaminated by bulges, bars, and 
non-circular motions. For a disk embedded in a dark 
halo, we expect Uz will increase for a given disk Edyn 
(Bottema 1993), hence the estimate is an upper limit. 

Combining Equation 1 with the expression for Sdyn, 
and normalizing k to an exponential density distribution 
and q for an oblateness of a typical disk with a 3.5 kpc 
scale- length (Equation 1 from DMS-II), we find the fol- 
lowing relation between the maximum rotation speed of 
the disk, disk central velocity dispersion, and disk oblate- 
ness: 
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Our expectations are, then, that for a maximum disk 
where V^'^x ~ 0.85Vinax, there should be a linear relation 
between Vmax and cr^'/jLo ''^ith a slope of ^0.43. 

However, in Figure la, the distribution of cr^'ljLo 
sus ^max for the DMS sample he near cr^'f^^o = (0.26 ± 
0.10)Kiax, a slope well below the expected value of 0.43 
for a maximum disk. The disks of our sample galax- 
ies would be maximal if hji/hz ~ 30. Based on the 
compilation from DMS-II (e.g., see Kregel et al. 2002), 
this oblateness is twice as flat as any observed edge-on 
disk. For realistic values of hn/hz, the disk-only rotation 
curves peak well below the measured T^nax- These galaxy 
disks are sub-maximal. 

To test this result, we plot the uncorrected central disk 
velocity dispersion (cr^og i?=o) filled gray points in 
Figure la. The corrections for the SVE projection are 
relatively small. Even without the correction disks still 
appear to lie in a parameter space that is marginally sub- 
maximal. This is an upper limit for any reasonable SVE 




V^ax (km s-') V,^ (km s"^) 

Figure 1. Stellar velocity dispersion versus maximum gas rotation speed: (a) it^'IJLq for face-on galaxies from the DiskMass survey (filled 

circles) and Bottema (1993; open circles); (b) '^'ji^^—ii^ for edge-on galaxies from Bottema (1993) and Kregel et al. (2005; open squares and 

MW for the Milky Way). Lines show the maximum rotational velocity allowed for a self- gravitating galactic disks of different oblateness 
(labeled). Maximum disks in the observed oblateness range would fall in the gray shaded region. 



with a < 1 and /3 within a factor of 2 of the epicychc 
value, and considering the above-mentioned impact of a 
dark halo on CTz. 
In Figure lb we plot tr^'*'^^^^ versus Vmax for the 

edge-on samples of Bottema (1993) and Kregel et al. 
(2005). As these authors state, for an exponential disk 
mass-density distribution of constant SVE with a — 0.6, 
'^R^R=hR ~ '^z^R=a- '^h^ similarity of the relation found 
for both edge-on and face-on samples demonstrates the 
accuracy of our basic assumption concerning the average 
SVE: If a were 65% larger on average (a = 1 instead of 
0.6), then the edge-on sample would be centered in the 
maximum-disk region of Figure 1. While this check is of 
primary importance for interpreting the maximality of 
edge-on disks, it also confirms our second-order correc- 
tions to (T^Q5 for face-on samples. A corollary, important 
below, is that the disk-oblateness measurements of these 
edge-on samples, used for mass-estimates of our face-on 
sample, are likely also an accurate application. 

4. WHY WE EXPECT ALL DISKS ARE 
SUBMAXIMAL 

We quantify disk- maximality as F^^^ = Klax/^max, 
using Equation 2 for V^f^ and Equation 1 from DMS-II 
for q. Figure 2 shows -F^iax of the DMS sample versus 
maximum rotation speed, if -band luminosity (Mk), and 
B — K color (velocity and luminosity are tightly corre- 
lated via itf)- We adjust q for the 3 galaxies in our 
sample earlier than Sb or later than Scd, as discussed 
in DMS-II. Our value of F^^l^ for UGC 463 agrees well 
within the errors to the detailed calculation presented 
in DMS-IV. For the DMS sample {F^lD = 0.47 ± 0.08 
(rms), equivalent to 22~^^% by mass for the disk within 
2.2hji, ignoring details of oblateness (the range is marked 
in Figure 2c). In contrast to Herrmann & Ciardullo 
(2009), we find little dependence of F^^^l^ on morphology 
(Figure 2c). There is also no significant trend of F^^^^ 
with disk central surface-brightness [{noiK)) = 17.9 mag 
arcsec"^, ranging from 16 to 20 mag arcsec"^; 18.1 mag 
arcsec"^ corresponds to a 'Freeman' disk for B — K = 3.5 



mag], bulge-to-disk ratio ((B/D) =0.1, ranging from to 
0.75), or bar classification. As Figure 2 shows, however, 
galaxies with higher rotation speed (see also Kranz et al. 
2003), greater near- infrared luminosity, and redder color 
are significantly more maximal. Adopting a regression 
model with intrinsic scatter (Akritas & Bershady 1996), 
we find i^^S = (0.24±0.08) + (0.26±0.08)(Kiax/200 km 
s-i) and F^ifk ^ (0.50 ±0.07) - (0.06 ± 0.02)(M;f-|-24), 
including bootstrap errors. The Vinax trend is consistent 
with disk mass- fractions doubling as (total mass)^/^. 

We have modeled the trend with color assuming galax- 
ies have a constant mean surface-brightness with color, 
and estimating changes in T with color based on all 
SPS models from Bell & de Jong (2001), Portinari et 
al. (2004) and Zibetti et al. (2009). Salient model vari- 
ations include treatment of chemical enrichment, star- 
formation history, and stellar evolutionary tracks (with 
varying amounts of TP-AGB stars). Different mass ze- 
ropoints (e.g., IMFs) are irrelevant here since we nor- 
malize the models at F^'fx = 0.525 ai B - K ^ 3.5. 
The curves in Figure 2c show the median and extrema 
of these model predictions for relative changes in T with 
color. The trends are upper limits since they assume dy- 
namical mass is equivalent to stellar mass, and discount 
contributions from disk gas. 

The dependencies of F^^^ on color, luminosity and ro- 
tation speed are such that even at extrema (maxima) in 
each quantity, the extrapolated value of F^^ < 0.75. 
For example, we have labeled the B — K colors for differ- 
ent Hubble and spectral types (Bershady 1995). Our 
sample spans essentially all of the color range. This 
means that even the earliest-type, fastest-rotating, red- 
dest disk should still be suhmaximal on average. In- 
deed, the most maximal disk from Herrmann & Ciardullo 
(M94; 2009) has a value of F^l^ - 0.7. Calculations for 
our sample assume k = 1.5, i.e., an exponential disk ver- 
tical density distribution. If, for example, disk vertical 
density distributions are isothermal, then our estimates 
of -Fniax decrease by ~15%; mass-fractions decrease by 
25%. 
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Figure 2. Ratio of disk rotation speed to total rotation speed (F^isi;) versus (a) total rotation speed, (b) K-band luminosity and (c) B — K 
color for DiskMass Survey galaxies. Maximum disks have Fjisij = 0.85 ±0.1 (shaded gray). Galaxies are coded by Hubble type (Sa, circles; 
Sb/Sbc, squares; Sc/Scd, triangles; Sd/Im, pluses) and bar classification (SA, blue/dark gray; SAB, green/light gray; SB, red/medium 
gray). Lines and labels are described in the text. 



These results have significant implications for the IMF, 
stellar evolution, cosmological accounting of baryons, 
and the formation of galaxy disks. Since a truncated IMF 
is already required to match T for a maximal disk (Bell & 
de Jong 2001), substantially sub-maximal disks require 
either uncomfortably top-heavy IMFs, or validates re- 
cent suggestions of a surfeit of luminous, low-mass stars 
(e.g., TP-AGB). A quantitative assessment of this sur- 
feit is forthcoming. Similarly, current accounting of the 
distribution of stellar mass in cosmological volumes is 
based on an T calibration that either assumes disks are 
maximal, or agrees with such a calibration (e.g., Li & 
White 2009, McGaugh et al. 2010). These stellar mass 
estimates should be lowered by at least a factor of 3 for 
disks, since represents the total dynamical mass, 

not just stellar mass. For example, the stellar mass of the 
disk of UGC 463 is only 60% of the total disk dynamical 



mass (DMS-IV). Unless the provocative claim for high 
T in elliptical cores (van Dokkum & Conroy 2011) is 
confirmed and shown to be more wide-spread, our down- 
ward revision of T is applicable to cosmological samples. 
Finally, galaxy formation models must reproduce accu- 
rate scaling-relations and their scatter (e.g., total mass 
to luminosity) in the context of submaximal disks. This 
submaximality must correlate with total mass, luminos- 
ity and stellar population, balancing trends of baryon 
loss from feedback with changes in angular momentum 
that, by our reckoning, on average must decrease three- 
fold from maximum-disk estimates. 
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